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FOREWORD 


This  report  has  been  written  as  Chapter  12 
of  a  forthcoming  book,  Engineering  Magnetohydro* 
dynamics.  References  to  other  chapters  refer  to 


other  chapters  in  that  book. 


12. 1  INTRODUCTION 

The  subject  matter  to  be  covered  in  the  present  chapter  will  consist  of 
Magnetohydrodynamic  flows  of  a  boundary  layer  character,  hi  conventional  fluid 
dynamics,  boundary  layer  theory  has  been  developed  to  a  high  degree  of  sophisti¬ 
cation.  No  comparable  extensive  treatment  will  be  attempted  here.  The  problems 
which  will  be  discussed  will  be  chosen  to  illustrate  those  new  features  introduced 
by  the  interaction  between  currents  flowing  in  the  conducting  fluid  or  plasma  and 
the  electromagnetic  field.  The  principle  interest  in  the  present  chapter  will  be 
in  flows  for  which  the  magnetic  Reynolds  number  is  vanishingly  small.  For 
completeness,  however,  the  boundeiry  layer  character  of  flows  with  R^  -*  «>  will 
be  treated  in  the  last  section. 

The  first  problem  considered  will  be  the  Rayleigh  problem  in  which  a.  magnetic 
field  is  applied  normal  to  the  surface  of  an  impulsively  moved  half  plane.  Its 
principle  value  will  be  in  that  a  solution  can  be  obtained  in  closed  form  so  that 
the  nature  of  Magnetohydrocfynamic  boimdary  layer  flows  can  be  inferred.  Follow¬ 
ing  this  introductory  treatment  the  boundary  layer  equations  are  obtained  from 
the  full  equations  by  means  of  the  well  known  boundary  layer  approximation. 

Based  on  these  equations,  incompressible  boundary  layer  flows  are  examined 
first.  Two  methods  of  solution  of  the  non-linear  partial  differential  equations 
are  developed.  For  flows  which  satisfy  certain  specified  conditions,  similarity 
exists  and  the  equations  can  be  reduced  to  an  ordinary  differential  equation  uiiich 
must  then  be  integrated  numerically.  For  more  realistic  boundary  conditions 
similarity  cannot  exist  so  that  in  this  instance  series  solutions  are  necessary. 

For  high  speed  flows  (hyperscmic)  compressibility  effects  must  be  considered,  so 
compressible  boundary  layers  are  considered  next.  In  this  case  the  boundary 
conditions  are  such  that  similarity  conditions  can  be  applied  and  the  solution  to 
the  problem  reduces  to  that  of  numerical  integration  of  two  coupled  ordinary 
differential  equations. 

To  ccmclude  this  chsq)ter  the  electric  current  and  induced  magnetic  field 
boundary  layers  which  can  form  on  surfaces  when  the  magnetic  Reynolds  number 
is  very  hi^  are  ccmsidered.  Althouf^  the  Magnetic  Reynolds  is  extremely  small 
in  almost  all  cases  of  interest  in  continuum  MagnetcAydrodynamics  the  above 
subject  is  of  interest  since  it  involves  application  of  old  boundary  layer  concepts 
to  a  new  physical  problem. 
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12.2  THE  RAYLEIGH  PROBLEM 

The  problem  which  will  be  studied  in  the  present  section  is  shown  in 
figure  12-1. 
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Figure  12-1.  Rayleigh  Problem  Configuration 


A  field  is  applied  perpendicular  to  the  plate  and  is  assumed  to  be  uniform  in 
space  and  constant  in  time.  The  doubly  infinite  plate  (a  non-conductor)  is  assumed 
to  move  in  x  direction  impulsively  at  t  =  0.  The  problem  is  to  determine  the 
induced  flow  and  magnetic  field  due  to  the  impulsive  motion  in  general. 

A  number  of  authors  have  investigated  the  general  problem^  and  have 
calculated  the  coupled  flow  and  induced  magnetic  fields  by  means  of  various 
approximations.  Of  principle  interest  in  the  present  chapter  will  be  flows  in 
which  R_  «  R„-  In  this  case  the  problem  reduces  to  one  in  which  the  flow 
disturbance  is  confined  to  a  thin  boundary  layer  region  while  the  induced  field 
varies  over  a  much  larger  region. 

The  equations  appropriate  to  the  present  problem  are  essentially  the  same 
as  those  for  the  transient  couette  flow,  except  that  now  the  second,  stationary, 
wall  is  infinitely-  far  away; 


P 


9u  _  ^ 

at  n 


®o  ay 


an 


(12.1) 


(12.2) 


The  boundary  and  initial  conditions  are 


t  =  0:  u  (y,  0)  =  0  b^^  (y,  0)  =  0 

t  >  0:  u  (0,  t)  =  u^  (0,  t)  =  0 
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where,  of  course,  u  and  must  remain  finite  as  y-*«.  The  requirement  that 
b  =  0  on  the  moving  non-conductive  wall  is  arrived  at  in  much  the  same  way 
as  was  done  for  Couette  flow,  see  Chapter  10  section  10.3. 

The  above  equations  can  be  most  readily  treated  with  the  aid  of  the  Laplace 
transform  defined  by 


L  |f  (y.  t)|  (y,  s)  =  /“e"®*  f  (y,  t) 


dt 


where  s  is  some  constant  parameter.  Taking  the  Laplace  transform  of  (12. 1) 
and  (12.2)  gives 
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(12.3) 
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o  dy  no  ^2 

(12.4) 

where  v  =  r]/p  =  kinematic  viscosity  and  the  boundary  conditions  are  now 

u  (o,  s)  =  bj^  (o,  s)  =  0 

The  solution  of  equations  (12.3)  and  (12.4)  is  readily  found  to  be 
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=  Ae-"^y  -H  Be-^y  (12.6) 

where  A  and  B  are  arbitrary  constants.  Terms  which  diverge  at  infinity  have 
been  omitted,  and  ±m  and  ±n  are  solutions  of 

-r*'  W  2  , 

-s)  (.r  -8)  - 

which  is  satisfied  if 
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The  constants  A  and  B  must  be  selected  to  satisfy  the  boundary  conditions 
imposed.  The  solution  then  assumes  the  following  form: 
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(12.9) 


Inversion  of  these  equations  is  quite  difficult  as  they  stand.  Since  our 
principle  interest  will  be  in  flows  in  which  —0,  and  this  implies 


or 


1 

op 


»i^  ,  a  solution  of  boundary  layer  character  can  be 


anticipated  and  the  problem  accordingly  simplified.  Keeping  finite,  and 
letting  0  while  holding  y/^  fixed  leads  to  the  following  results: 
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Accordingly,  when  =  0  at  the  moving  wall  it  is  zero  throughout  the  velocity 
boundary  layer  region.  Such  a  result  could  have  been  anticipated.  If  on,  the 
other  hand,  other  boundary  conditions  had  been  selected  for  b^  other  results 
would  have  been  obtained.  If  the  problem  geometry  were  such  that  for  long 
times  ,  b  (0,  t)  =  constant  rather  than  zero  then  the  boundary  condition  would  be 
that  b^  at  the  boundary  be  some  function  of  time  depending  on  the  total  current 
flow.  Such  a  boundary  condition  could  not  be  conveniently  handled.  Finally,  if 
the  moving  wall  is  assumed  to  be  a  perfect  conductor  the  electric  field  in  the 
fluid  adjacent  to  the  wall  must  be  zero.  Then  from  equations  (10.1)  and  (10.5a) 
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we  have 


r  1 


or 


"87  (o.  t)  =  -  (T  Uoo  B^ 


(12.11) 

which  is  the  required  boundary  condition  for  this  case.  As  has  been  shown  by 
Ludford^^^  the  solution  when  equation  (12. 11)  is  used  as  the  boundary  condition 
leads  to  a  b  which  is  constant  throughout  the  boundary  layer  region. 

A 

Returning  now  to  the  problem  originally  posed,  the  inverse  of  equation  (12. 10) 
is  readily  found  to  be 
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This  result  is  identical  to  that  of  Rossow'  ’  although  the  induced  magnetic  field 
was  not  considered  in  that  reference.  The  reason  for  the  agreement  is  obvious 
as  bx  =  0  throu^out  the  boundary  layer  region.  However,  it  should  be  noted 
that  bj^  will  become  finite  above  the  boundary  layer  and  approach  some  constant 
value  at  «>.  The  value  of  formulating  the  complete  problem  before  making  the 
boundary  layer  approximation  will  be  seen  shortly  when  the  shear  stress  at  the 
wall  is  calculated. 

First,  however,  consider  the  solution  for  u  found  above, 
to  go  to  zero  equation  (12.12)  reduces  to  the  classical  Rayleigh  problem 


If  is  allowed 
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which  is 


In  this  case  the  velocity  profile  depends  on  only  one  variable 

2\/1T 

essentially  the  similarity  variable  to  be  discussed  later.  In  the  presence  of 
magnetohydrodynamic  effects  such  similitude  no  longer  exists  and  there  is  no  one 
variable  on  which  depends,  hi  order  to  illustrate  the  flow  pattern,  calcu- 

lations  based  on  equation  (12.12)  have  been  carried  out  and  are  presented  in 
Figure  12-2. 
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Figure  12-2.  Velocity  Profiles  in  the  Rayleigh  Problem 

It  can  be  seen  that  application  of  the  magnetic  field  increases  the  time  required 
for  the  flow  velocity  above  the  plate  to  reach  any  specific  value. 

If  desired,  the  above  solution  can  be  interpreted  as  a  boundary  layer  flow 
with  t  -*  x/uoo.  In  this  case  we  observe  that  as  the  field  strength  is  increased 
the  velocity  profile  becomes  fuller.  This  result  is  analagous  to  that  found  for  the 
Hartmann  flow  earlier.  Also,  the  above  results  suggest  that  a  natural  variable 
to  formulate  the  boundary  layer  problem  with  should  be  t)  =  —  just  as  in  the 

_  /  IfV 

in  the  case  of  conventional  incompressible  boundary  layers. 
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Finally,  the  force  on  the  moving  plate  will  be  considered, 
is  given  simply  as 


(^)co  =  ^ 


The  shear  stress 


(12.14) 


Now  from  equation  (12.8)  the  transform  of  may  be  derived 
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The  inverse  of  this  expression  gives  the  following  result 


2 

where  ctB^ 

y  =  - ^  V. - 

p  (1+v^ ) 

which  agrees  with  Rossows^^^  result  when  P  =  0.  The  above  equation  is  valid 
for  arbitrary  values  of  and  is  thus  a  much  more  general  result. 

hi  order  to  illustrate  the  influence  of  on  the  solution,  equation  (12. 16) 
is  plotted  in  figure  12-3. 


A  number  of  valuable  observations  can  be  made  based  on  these  results. 
First,  it  is  clearly  seen  that  as  a  increases  the  wall  shear  becomes  larger 
for  any  given  time.  This  increased  shear  arises  as  a  result  of  the  enhanced 
Lorentz  force  on  the  conducting  fluid  near  the  moving  wall.  At  very  long  times 
all  of  the  curves  of  Figure  12-2  reach  asymptotic  values  given  by 


lim  t 


«  (Tw') 


Now  in  the  non-magnetohydrodynamic  Rayliegh  problem  equation  (12. 13)  shows 
that  u-*«ast-*ooso  that  -*  0,  and  it  is  of  interest  to  inquire  as  to  why 
in  the  present  problem  ^  0  in  the  same  limit.  If  we  return  to  the  original 

o 

equations  of  the  problem,  equations  (12.1)  and  (12.2),  and  assume  =  0  the 
following  solution  is  obtained  for  u  and  b^. 


u 

Uoo 

bx= 

The  steady  state  long  time  solution  then  has  an  exponential  velocity  profile  rather 
than  the  constant  one  when  say  0  =  0.  This  shape  of  the  profile  arises  because  of 
the  existence  of  the  currents  in  the  fluid  and  the  attendent  Lorentz  force.  The 
form  of  the  current  distribution  is  also  exponential  decreasing  from  some  finite 
value  at  the  wall  to  zero  at  infinity. 
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From  the  point  of  view  of  the  physical  problem  the  wall  shear  is  balanced 
by  the  reaction  force  on  the  source  of  the  magnetic  field.  When  there  is  no 
magnetic  field  there  is  nothing  to  balance  this  shear  and  so  it  must  be  zero  in 
the  steady  state. 

12.3  FORMULATION  OF  BOUNDARY  LAYER  EQUATIONS 

The  boundary  layer  approximation  in  fluid  mechanics  has  been  discussed  by 
many  authors  from  the  time  of  Prandtl;  there  is,  therefore,  no  justification  in 
repeating  these  arguments.  What  will  be  discussed  in  the  present  section  will 
be  those  new  terms  which  appear  in  the  equations  due  to  Magnetohydrodynamic 
effects,  and  their  simplification  via  the  boundary  layer  approximation.  Attention 
for  the  present  will  be  restricted  to  flows  with  R^  ~  0  so  that  the  magnetic 
field  can  be  taken  as  the  applied  field  (i.e.  section  10.8  of  chapter  10).  The 
question  of  the  boundary  layer  approximation  when  R^  —  will  be  discussed 
later. 

The  general  problem  which  will  be  considered  is  shown  in  Figure  12-4. 

The  magnetic  field  vector  is  assumed  to  lie  in  the  xy  plane  but  may  have  x  and 
y  components.  The  Lorentz  force  for  such  a  two  dimensional  problem  is  given 
by 

F  =  2  X  =  ^[e  +  q  X  b]  X  B  (12. 17) 


B 


Figure  12-4.  Magnetohydrodynamic  Boundary  Layer  Problem 
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where  the  current  is  derived  from  the  Ohm's  law  since  induced  magnetic  fields 
have  been  assumed  negligible.  Now  since  only  steady  flow  problems  will  be 
considered  it  will  be  valid  to  take  E  =  0  or  if  desired  some  constant.  It  may 
be  recalled  that  in  the  unsteady  problem  (i.e.  transient  Couette  flow)  it  is  not  in 
general  valid  to  assume  E  =  0.  For  the  present  case,  equation  (12.17)  can  be 
written  as 

(12. 18) 


F^=-<T(uBy 


-  By) 


Fy= 


vBx^) 


(12.19) 


For  purposes  of  order  of  magnitude  arguments  one  can  assume  that  B  =  0(B  ). 

V  y  ^ 

Then,  since  —  =  0  (6)  where  6  is  the  boundary  layer  thickness  the  first  of  the 

above  expressions  can  be  written  as 


ff  B  ^  u 
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(12.20) 


where  By  can  be  a  function  of  x  and  y  in  general. 


Li  a  boundary  layer  calculation  it  will  be  exceedingly  inconvenient  to  keep 


By  a  function  of  y. 


Accordingly,  it  is  of  interest  to  study  By  (x,  y)  somewhat 


more  carefully, 
requires 


For  a  two  dimensional  field  one  of  Maxwell's  equations  (V*  B  =  0) 
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Now  with  B^  =  0(B^^)  and  •^  =  6  (1)  while  •^  =  0  (6"^)  we  have 


9y 


9B, 


=  0(1) 


or 


e  (6) 


Thus,  the  change  in  By  across  the  boundary  layer  will  be  the  order  of  6  and  can 
be  neglected. 

Making  use  of  equation  (12.20),  assuming  p  =  p(x)  and  making  the  usual 
boundary  layer  approximations  the  x  component  of  the  momentum  equation  becomes: 
becomes: 


P  (u 
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Before  proceeding  further  the  assumption  that  p  =  p(x)  must  be  examined 

further.  Again  assuming  B  =  0  (B  )  equation  (12. 19)  reduces  to 

y  " 
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Now  the  y  component  of  the  momentum  equation  can  be  written  in  the  following 
form 
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where 


Q  = 
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and  the  velocities  have  been  made  dimensionless  by  some  reference  velocity  (Uoo), 
the  distances  x  and  y  by  some  reference  length  (L),  the  pressure  by  PqUoo  ,  the 
density  p  and  viscosity  tj  by  reference  values,  and  the  magnetic  field  by  some 
reference  field  (Bj^).  For  magnetohydrodynamic  boundary  layer  problems  of 
interest  Q  =  0(1).  Then  an  order  of  magnitude  analysis  of  equation  (12.23)  with 
u  =  G(l).  V  =  0  (6).  ^  =  0(1),  |r  =  0  (|).  and  Rg  =  0  (-I)  gives 

=  0(1) 

or 

Ap  =  0(6)  ~  0 


Thus,  the  pressure  change  across  the  boimdary  layer  is  the  order  of  6  (a  small 
quantity)  and  can  be  neglected.  It  should  be  noted,  however,  that  in  the  absence 
of  magnetic  forces  A  p  =  0  (6  )  so  that  it  is  less  valid  to  assume  p  =  p  (x)  in 
the  present  case  than  it  had  been  in  the  absence  of  a  magnetic  field. 

Next,  the  Joule  heating  term  which  will  enter  into  the  energy  equation  must 
be  considered. 

.2 

H  =  ^  =  0^  (3  X  B)  •  (q  X  B)  (12.24) 
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as  before,  the  last  two  terms  can  be  neglected  since  v/u  =  0(6).  Thus, 
H  s  (T  u^  By^ 
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and  the  energy  equation  can  be  written  as  follows 
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Finally,  the  mass  continuity  equation  must  be  added  to  complete  the  system 
of  equations: 


h  (pu)  +  ^  (pv)  =  0 


(12.26) 


Before  proceeding  to  the  solution  of  specific  problems  it  will  be  of  interest 
to  study  the  general  question  of  similar  solutions  for  incompressible,  constant 
property,  boundary  layer  flows.  Consider  equation  (12.21)  evaluated  in  the  inviscid 
free  stream  flow 
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Combining  this  with  equation  (12.21)  leads  to  the  following  relation. 
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Introducing  u  =  -|^  ,  and  v  =  -  -|^  satisfies  continuity  and  equation  (12.28) 
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T  yx 


X  yy 


+  V  ip. 


yyy 


(12.29) 


The  problem  now  is  to  determine  whether  or  not  the  above  can  be  reduced  to  an 
ordinary  differential  equation  by  a  proper  choice  of  variables.  Let  us  assume 


Ip  (x,y)  =  y'uoo  vx  f  (n) 

where 

V  (x,y)  = 

and  where  it  should  be  noted  that  uw  and  B^  are  known  functions  of  x  as  yet 
unspecified.  Substituting  the  above  equation  (12.29)  becomes 
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2  ^  (f)^  -  ff" 

2 

/I  A.ft\  =  4-  ^22  fill 

p  -  2  ^  ^  “ax~  X  ^ 

Accordingly,  similarity  can  be  achieved  if 

dUoo  ^  Uoo 

dx  X 


and 


(12.30) 


so  that  Uoo  and  By  must  be  of  the  following  form  to  permit  a  similar  solution 

Uoo  =  Cj^  X™ 

m-1 

B  =  Co  X  ^ 

y  2 

In  fluid  mechanics  the  x™  variation  of  Uoo  corresponds  to  the  wedge  flow  solutions. 
In  magnetohydrodynamic  problems  such  a  simple  interpretation  is  not  possible 
since  the  applied  magnetic  field  interacts  with  the  inviscid  free  stream  flow. 

We  will  return  to  this  point  in  the  next  section. 

A  final  point  in  regard  to  the  above  treatment  should  be  made,  hi  looking 
for  similar  solutions  only  the  incompressible  boundary  layer  was  studied.  For 
a  compressible  flow  the  situation  is  much  more  complex.  Now,  for  a  similar 
solution  to  be  found  for  the  compressible  boundary  layer  problem  it  is  safe  to 
say  that  a  minimum  condition  would  be  that  Uoo  and  satisfy  the  conditions 
already  established  for  the  incompressible  case,  hi  addition,  many  other 
assumptions  and  requirements  will  be  necessary  to  obtain  similarity  in  the 
compressible  case.  Some  of  the  new  phenomena  associated  with  compressible 
magnetohydrodynamic  boundary  layer  flows  will  be  treated  in  section  12.5. 

The  following  section  will  be  restricted  to  incompressible  boundary  layers. 

12.4  INCOMPRESSIBLE  BOUNDARY  LAYERS 


The  distinguishing  feature  of  an  incompressible  magnetohydrodynamic 
boundary  layer  is  that  the  inviscid  flow  external  to  the  boundary  layer  is  also 
a  conductor.  Then,  since  it  is  not  possible  to  restrict  the  magnetic  field  to 
the  boimdary  layer  region  alcHie,  one  has  no  choice  but  to  incorporate  the 
results  of  the  inviscid  magnetohydrodynamic  analysis  into  the  boundary  layer 
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investigation.  On  the  other  hand,  this  may  not  be  necessary  in  compressible 
flows  in  which  the  conductivity  can  vary,  and  for  certain  types  of  flows  may  be 
assumed  zero  external  to  the  boundary  layer.  Application  of  the  similarity 
solution  to  incompressible  boundary  layer  although  feasible,  in  principle,  will  not 
be  carried  out  due  to  the  difficulty  of  interpreting  the  resulting  inviscid  flow. 
Instead,  two  simple  incompressible  boundary  layer  problems  will  be  analyzed  by 
the  series  expansion  technique 

Boundary  Layer  with  Uniform  Free  Stream 

The  first  problem  that  will  be  considered  will  take  as  a  model  for  the 
inyiscid  flow  the  flow  through  a  parallel  walled  two  dimensional  channel  and  will 
consider  the  flat  plate  situated  somewhere  in  the  channel  some  distance  from 
either  wall,  hi  this  case,  the  free  stream  velocity  will  be  a  constant  and  the 
pressure  gradient  will  be  given  by  dp/dx  =  -ouoo  .  The  applied  magnetic 
field  will  be  assumed  constant  so  that  the  pressure  will  vary  linearly  along  the 
channel.  The  equation  to  be  solved  is  equation  (12.29)  which  simplifies  to  the 
following.  2 

lb  lb  -  lb  lb  +  - ^  {ip  -  Uoo)  =  u  ip  (12.31) 

^y  ^yx  ^x  ^yy  p  '^y  ’  ^yyy  '  ' 


where  the  boundary  conditions  are 

at  y  =  0 
at  y  =  «,  X  =  0 

Since  the  pressure  gradient  in  the  present  problem  is  constant  and  favorable  it 
will  be  adequate  to  use  a  classical  Blasius  series  expansion  to  effect  a  solution. 


ip  =  ip  =0 

ry 

</'y  =  Uoo 


Assume 


_ _  r  ctB  2  /oB  M2  1 

ip  (x,y)  =yu«,yx  [fo  +  x  f^  +  ^  xj  fg  +  . . .  .J 


(12.32) 


where  it  is  assumed  that 


a  small  quantity  and  the  f's  are  functions 


of  Tj  =  •  Introduction  of  the  above  assumption  for  ip  into  equation  (12.31) 

leads  to  the  following  infinite  set  of  ordinary  differential  equations. 

f  •"  +  f  f '*  =  0, 

o  o  o 

with  the  following  boundary  conditions: 


(12.33) 


fo  (o)  =  ^o'  ® 


V 
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(12.34) 


with 


h  2  ^1  2  *0  h  ^ 


(o)  =  (o)  =  0  (oo)  =  0 


f  in  4-  i.  f  f  M  _  9f  I  f  I  4-  ^  f  II  f  -  It  i\^  _  ^  f  f  II  +  f  I 

h  ^2  ^0^2  ^2  ^  2  h  '^1  )  2  hh  4 


(12.35) 


with 


f„  (o)  =  f  ’  (o)  =  0 


f  '  (oo)  =  0 


etc. 

The  first  of  these  is  the  well  known  Blasius  equation  for  which  the  solution 
has  been  tabulated  for  the  boundary  conditions  shown.  All  subsequent  equations 
are  linear,  but  since  they  depend  on  the  preceeding  solutions  numerical  integration 
is  necessary.  The  first  two  equations  beyond  the  Blasius  have  been  solved  by 
Rossow^^\  The  resulting  velocity  profiles  are  shown  in  figure  12-5. 


Figure  12-5.  Incompressible  Boundary  Layer  for  Uniform  Applied 
Magnetic  Field  and  Free  Stream  Velocity 

Again  these  results  are  in  agreement  with  the  Hartmann  flow  treated  in 
Chapter  9  and  the  Rayliegh  problem  solution  discussed  earlier  in  the  present 
chapter.  In  other  words,  the  presence  of  the  applied  magnetic  field  tends  to 
make  the  velocity  profile  fuller.  If,  on  the  other  hand,  the  free  stream  velocity 
had  been  decreasing  with  x  the  pressure  gradient  would  have  become  unfavorable, 
the  boundary  layer  profile  would  have  become  less  full,  and  the  possibility  of 
flow  separation  would  exist.  Such  a  case  will  be  considered  next. 
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Boundary  Layer  Subject  to  Adverse  Pressure  Gradient 

Consider,  as  an  appropriate  inviscid  flow  the  flow  through  a  parallel  walled 
channel  when  the  applied  magnetic  field  is  created  by  a  current-carrying  wire 
alligned  perpendicular  to  the  flow  direction  and  imbedded  in  the  lower  wall.  As 
shown  in  Figure  12-6  the  boundary  layer  development  along  the  lower  wall  will 
be  considered. 


Figure  12-6.  Boundary  Layer  Development  with  Magnetic  Field 
Generated  by  Current-Carrying  Wire  in  Lower  Wall 

One  must,  of  course,  assume  that  the  boundary  layer  ahead  of  x  =  0  can  be 

removed,  perhaps  by  a  bleed  port  shown  schematically  as  the  dotted  region. 

The  equation  and  boundary  conditions  for  ip  have  already  been  presented, 

equation  (12.31),  and  will  also  be  used  in  the  present  problem.  One  should  be 

aware,  however,  that  requiring  =  u*  at  y  =  «>  means  that  the  inviscid  flow 

shear  (see  Figure  10-33)  at  the  lower  wall  is  being  neglected.  In  general,  this 

shear  can  be  neglected  when  the  inviscid  free  stream  flow  has  been  linearized. 


As  noted  earlier,  the  existence  of  a  non-uniform  free  stream  velocity,  and 
the  possibility  of  separation,  requires  the  use  of  a  series  solution  which  is  more 
sophisticated.  A  procedure  frequently  used  for  conventional  boundary  layers 
and  ideally  suited  to  the  present  problem  is  due  to  Gortler.  Its  extension  to 
magnetohydrodynamics  and  its  application  to  the  present  problem  will  be  carried 
out  now. 

If  the  following  dimensionless  variables  are  defined 


=  2SS. 


=  Z. 


Uoo 


Uoo 

U. 


iPyfR^ 

y  iTT 
*^o  1 


B 

B  =  ^ 

y  ®o 
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and 


%  =  ^  yo  " 


Q  = 


^  y 
o  0 


P  U 


i 


where  =  u  (o,y),  assumed  constant,  y^  is  some  reference  length,  and 
some  reference  magnetic  field,  and  the  following  new  independent  and  dependent 
variables  are  defined 


«=  r 

o 


Uoo  dx 


=  Uoo  y 

yfu 


F  =  ^  /yfu 

Then  equation  (12.31)  now  becomes 


F  +  F  F  + 
TjT)T]  rrn 


'ii^l 
u„2  J 


(1  - 


=  2| 


B. 


~  ^  ~  2 
Uoo 


(F.  -  1) 


] 


with  the  following  boundary  conditions 

at  T)  =  0 
at  Tj  =  oo 


F  =  F„  =  0 
V 


(12.36) 


F„  =  1 

V 

The  problem  can  be  solved  in  general  by  assuming  the  following  expansions; 
F  a.V)  =  F^  (T,)  +  ^  Fj  (Tj)  +  Fg  (Tj)  +  ... 

2 


Uoo 


=  ^(1)  =  0+  ^0,  +  i^0^  + 

dx  0  12 


•  •  •  • 


B  ^ 

=  g  (0  =  go  +  I  gj  +  g2  - 

Uoo 

Substituting  these  into  equation  (12.36)  one  obtains,  as  before,  an  infinite  set  of 
equations.  The  first  of  these  is  the  Blasius  equation 

V  ^o  V  "  ®  (12.37) 

Fq  (o)  =  F^'  (o)  =  0  F^'  (CO)  =  1 

The  second  and  all  later  equations  are  linear,  and  are  given  by  the  following 
recursive  system 


-17- 


Fy  +  F^Fj^’ '  -  2kFJjFJ^  +  (2k  +  1)  F|j'  Fj^ 


(12.38) 


where 


^  ^-1  ®  ^  1.  2,  3 . 


,  k-1 

k-j 

«k-i 

2 

i=l 

^j-i  n  '■k-i-i 

k-1  k-1 

-  ^Ic-i  -  2  jii 

(k 

-  «  Fj 

k-1 

+  _  S  (  1  +  2j  )  F.  F”_.  +  2g^_^  (1  -  f;^) 
k-1 

j=l  ^3-1  k-j 

and 

Fk  =  Fj^  =  0  at  Tj  =  0 
F’ 

k  =  0  at  T]  =  «. 

Application  of  the  above  equations  to  the  solution  of  the  specific  problem  of 
Figure  12-6  has  been  carried  out  by  Sherman Numerical  results  were  obtained 
for  an  inviscid  flow  in  wliich  L/y^  =  1.6  and  Q  =  0.1.  Curves  showing  the  wall 
skin  friction,  boundary  layer  thickness,  and  velocity  profiles  are  shown  in 
Figures  12-7  and  12-8. 


Figure  12-7.  Wall  Skin  Friction  Versus  Distance  From  Leading  Edge 

for  Q  =  0. 1  and  L/yo  =  1.6.  Nonmagnetic  case - ; 

Magnetic  Case,  Series  Solution - ;  Magnetic  Case, 

Integral  Approximation - • - . 
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Figure  12-8.  Boundairy  Layer  Velocity  Profiles  at  Several  Positions 
Downstream  of  Leading  Edge  for  Q  =  0. 1  and  L/y^  =  1.6 

In  order  to  interpret  these  results  properly  it  will  be  of  value  to  review 
the  physical  phenomena  occuring  in  the  flow.  First,  the  applied  magnetic  field 
acts  on  the  inviscid  flow  in  such  a  way  as  to  cause  the  boundary  layer  free  stream 
velocity  to  decrease  rapidly  in  the  vicinity  of  the  wire.  This  tends  to  make  the 
velocity  profile  less  full,  and  thereby  reduce  the  wall  skin  friction.  In  addition, 
the  magnetic  field  within  the  boundary  layer  creates  a  Lorentz  force  which  tends 
to  retard  the  flow.  This  effect  then  also  tends  to  reduce  the  skin  friction,  and 
both  combined  may  retard  the  boimdary  layer  sufficiently  to  cause  separation. 

That  such  flow  separation  is  indeed  a  practical  possibility  is  shown  in  Figures 
12-7  where  t  =  0  at  x  =  1.22  L,  just  beyond  the  wire.  For  the  present  strongly 
retarded  flow  the  classical  Karman-Pohlhausen  integral  approach  is  seen  to  be  a 
poor  approximation.  The  details  of  the  velocity  profiles  are  shown  in  Figure  12-8 
where  it  must  be  kept  in  mind  that  uw  is  decreasing  as  x  increases  from  zero. 

12.5  COMPRESSIBLE  BOUNDARY  LAYERS 

Although  some  of  the  requirements  for  similar  magnetohydrodynamic  boundary 
layers  had  been  identified  earlier  such  solutions  were  not  sought  in  the  incom¬ 
pressible  case  due  to  the  difficulty  in  interpreting  the  results  in  terms  of  a 
practical  problem.  When  the  flow  is  compressible,  however,  such  difficulties 
may  not  exist.  For  example,  when  considering  the  hypersonic  flow  over  a  flat 
plate  in  which  the  free  stream  is  at  a  low  temperature  compared  to  the  hi^ 
temperatures  in  the  boundary  layer,  one  may  assume  (Too  =  0  so  that  the  applied 
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magnetic  field  does  not  disturb  the  inviscid  flow  and  Uoo  =  constant.  Also,  there 
are  some  practically  important  problems  for  which  the  interaction  between  the 
magnetic  field  and  the  inviscid  flow  have  been  calculated  and  have  shown  Uoo  to 
have  the  required  form  for  similarity. 

In  the  cases  just  cited  similar  solutions  have  practical  significance,  and  in 
fact  are  absolutely  essential  for  progress  to  be  made  in  solving  these  complex 
problems.  In  the  present  section  two  compressible  boundary  layer  flows  will  be 
treated  in  some  detail,  and  a  third  will  be  discussed. 

First,  consider  the  hypersonic  flow  over  a  semi-infinite  flat  plate  when  the 
temperature  within  the  boundary  layer  is  high  enou^  to  ionize  the  gas.  In  order 
to  obtain  a  similar  solution,  the  free  stream  velocity  is  assumed  constant  (low 
free  stream  temperature  so  that  0=0),  and  the  applied  magnetic  field  is 

assumed  to  vary  inversely  as  the  square  root  of  x.  In  addition,  the  wall  tempera¬ 
ture  is  assumed  constant,  and  the  gas  is  considered  to  be  in  thermodynamic 
equilibrium.  The  geometry  and  coordinates  of  Figure  12-4  will  be  used. 

Neglecting  heat  flux  due  to  diffusion  of  species  (current  does  not  flow  across 
a  temperature  gradient),  the  Hall  effect,  and  induced  magnetic  fields,  and  putting 
^  =  0  the  boimdary  layer  equations  are  given  by  Equations  (12.21),  (12.25),  and 
(12.26).  The  magnetic  field  can  be  expressed  as 


The  reduction  of  these  equations  to  the  ordinary  differential  equations  corresponding 

(41 

to  similarity  is  accomplished  by  the  Crocco  transformation'  '.  The  procedure  will 
be  sketched  briefly. 

If  the  independent  variables  of  the  problem  are  changed  from  x  and  y  to  x  and 
9u 

u  and  T  =  M  ^  the  transformation  formulas  can  be  written  as  follows: 

(4 -(4 -(site). 

recognizing  that  we  can  consider  y  =  y  (x,  u)  the  first  of  these  leads  to 
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or 


and  the  second  yields 


or 


yu  = 


1 

T 


Then  the  transformation  equations  can  be  rewritten  to  yield 


Using  these  relations  the  transformed  boundary  layer  equations  become 

9  (Pv)  =  9  (Pu)  _  Ji  d  (pu) 

9u  9u  ~  T  9x 


pv  =  (pu)  y^  + 


dr 

9u 


_u 

T 


(12.39) 

(12.40) 


pu 


9h 

9x 


9u 


+ 


(12.41) 


where  the  second  term  in  the  last  of  these  has  been  simplified  with  the  use  of 
equation  (12.40).  Eliminating  pv  from  equation  (12.40)  by  use  of  equation  (12.39) 
gives 


pu 


+ 


+ 


2  ^  ,gqu 

o  S  du  '  T  ' 


=  0 


or 


(p-“)  * 

X 


Next  assume  that 

T  (X,  u)  = 


Gjul 


du 


ZIULx  = 


(12.42) 


and  h  (x,  u)  =  h  (u) 
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so  that  equations  (12.41)  and  (12.42)  become 


or 


and 


PUT]  1  1  d^G  ^o  ^0  d  /  g7?u  \  _ 

2°  ,/T  ^  du^  ■  VT  G  ;  ■  " 

_  d^G  _  PUT)  p  2-  d  /  CTT)u\  „ 

°  ^  ■  2  ^  ^®o  ^o  ^\—)  <12.43) 


G_ 

dG 

L  u 

0  0 

T)X 

du 

du 

X 

du 

2„  2 
dh\  ^  ®o 

^0 

T)X 

du 

w 

du  y  X 

.  + 

T)X 

_d 

Ur  du/ 

r' 

du 

du 

du 

0 

2  - 
uLqT) 

F—  +  ul 

G 

Ldu  “j 

These  are  then  two  ordinary  differential  equations  for  G  and  h  as  functions  of  u. 
The  numerical  procedure  for  their  solution  is  given  in  detail  by  Bush<®^  where 
the  method  of  expressing  p,  t),  cr,  and  as  functions  of  h  for  high  temperature 
air  are  also  described.  Some  of  the  results  of  calculations  for  the  constant  wall 
temperature  case  are  shown  in  Figure  12-9. 


Q 


Figure  12-9.  Shear  Stress  and  Heat  Flux  for  Hypersonic 

Boundary  Layer 
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The  conditions  chosen  for  the  above  case  were  Too  =  222®K,  p  =  10  atm, 

Moo  =  25,  and  a  wall  temperature  of  ~2000‘’K.  The  subscript  NM  stands  for 

—  2 
non-magnetic,  and  Q  is  the  interaction  parameter  proportional  to  . 

The  above  results  bear  a  striking  resemblance  to  the  real  gas  Couette  flow 
results  described  in  Chapter  10.  Again,  a  hysterisis  effect  exists  due  to  the 
form  of  the  a  versus  T  curve  and  the  dotted  portion  of  the  curves  are  unstable. 
The  principle  new  information  here  is  that  in  allowing  the  boundary  layer  to  grow 
in  height  (in  the  couette  flow  6  is  fixed)  a  reduction  in  heat  flux  of  80%  is  possible 
where  little  if  any  heat  flux  reduction  had  been  predicted  earlier.  The  shear 
stress  does,  however,  behave  as  predicted. 

Another  compressible  boundary  layer  flow  of  practical  interest  is  one  which 
grows  along  the  electrode  surface  of  a  crossed  field  MHD  channel.  The  principle 
new  feature  here  is  that  a  current  flows  normal  to  the  boundary  layer  surface 
(see  Figure  12-10)  so  that  the  contribution  to  the  heat  flux  due  to  the  diffusion 
of  electrons  in  a  temperature  gradient  cannot  be  neglected. 


Figure  12-10.  Electrode  Boundary  Layer  Configuration 


Since  the  boundary  layer  imder  consideration  will  be  developing  within  a  finite 
width  channel  the  usual  assumption  that  the  boundary  layer  thickness  is  not  large 
enough  to  disturb  the  inviscid  flow  must  be  made.  In  addition,  it  is  necessary 
to  assume  that  the  electrical  resistance  of  the  boundary  layer  is  small  compared 
to  the  resistance  of  the  inviscid  flow  so  that  the  overall  current  flow  is  determined 
external  to  the  boundary  layer. 

The  momentum  equation  is  given  by  equation  (12.21)  except  that  now  the 
Lorentz  force  is  written  as  j  B^,  since  (j)  is  now  some  known  function  of  x. 


pu 


‘•'’Si* 


dp  _  9 

dx  ~  ay 


ay) 


iB, 


(12.45) 
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The  energy  equation  written  in  terms  of  temperature  (h  =  c  T)  is  given  by 

P  2  / 

equation  (12.25)  except  that  here  the  Joule  heating  is  written  as  j  /o  again  due 
to  the  fact  that  j  =  j(x)  is  a  given  quantity.  The  continuity  equation  is,  of  course, 
equation  (12.26). 

Since  (j)  is  the  same  within  the  boundary  layer  as  in  the  free  stream  the 
momentum  equation  in  the  free  stream  is 

^  ^  “  J  ®y  <12-46) 

and  combining  this  with  equation  (12.45)  gives 


pu 


ax 


+ 


a  /_  9u  V 
ay  ay^ 


■<"  Poo  Uoo 


dUoo 

dx 


(12.47) 


The  energy  equation  evaluated  in  the  free  stream  yields 


poo  Uoo 


C 


P 


dToo 

dx 


Uoo 


^  + ii 

dx  (Too 


(12.48) 


which  when  combined  with  equation  (12. 25)  yields  the  following  relation 


(12.49) 


It  is  interesting  to  note  at  this  point  that  the  momentum  equation  for  the  present 
case  is  independent  of  B^.  The  principle  Magnetohydrodynamic  effects  appear  in 
the  energy  equation. 

For  convenience  in  the  boimdary  layer  analysis  n  and  K  will  be  assumed 

proportional  to  T,  despite  the  fact  that  for  constant  mean  free  path,  kinetic 
1/2 

theory  leads  to  a  T  dependence,  hi  general  (see  chapter  5)  the  electrical 
conductivity  can  be  determined  as  a  function  of  pressure  and  temperature.  This 
time,  for  ccmvenience  in  both  the  free  stream  and  boundary  layer  analysis  the 
pressure  dependence  will  be  ignored.  Finally,  the  gas  will  be  assumed  perfect 
so  that 


p  =  pRT 


(12.51) 


Now,  the  boundary  conditions  which  are  needed  are 


Too  =  Too  (X)  Uoo  =  Uoo  (x) 

hi  addition,  expressions  are  needed  for  p  (x)  as  well  as  j  (x).  All  four  of  these 
relations  can  be  obtained  from  the  solution  of  the  inviscid  problem.  For  the 
inviscid  flow,  it  will  be  assumed  that  Too  and  Eoo  are  constant,  the  former  also 
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leading  to  a  constant  (Too  if  the  pressure  dependence  of  ct  is  neglected.  Within 
these  assumptions  the  solution,  based  on  the  methods  of  chapter  11,  is 

Uoo  = 


P  = 


CToe  (RToo  Eoo) 

2 

a 


(5n-l)  l-5n 

2  ^ 
n 


(12.52) 


]  =. 


RTo 


2 

0!  n 


(5n-l)  X 


-2n 


where  one  must  have  n  >  1/5. 

As  noted  earlier,  Uoo  in  the  above  form  may  make  a  similarity  solution 
possible.  The  feasibility  of  reducing  the  equations  to  similar  form  will  be 
taken  up  next.  Define  the  following  independent  variables 

I  =  ^  ^  ^  (12.53) 

o  Pq 


and 


(12.54) 


where  (  )jj  denotes  some  convenient  reference  x  position.  Next,  define  a  stream 
function  in  order  that  mass  continuity  be  satisfied 


ij)  = 

y  pQ 


X  Po 


and  then  redefine  the  stream  function  to  be 


The  momentuin  equation  then  becomes 

f  T I  f  -j.  f  £  f  I  4.  ^  ^ 

Uoo  dx 

If  a  dimensionless  temperature  is  defined  as  0  =  the  energy  equation 

becomes  ^ 


=  0 


(12.55) 


1  d^e  f  ^ 


(y-1)  Moo^  (f’)^ 


CpPoUoo  p  ae 


-25- 


(12. 56) 


2u„  « 


Cp  Too  poo  Uoo 


—  0  (  ^  -  f  ) 

p  '  CJ  ' 


Using  the  results  of  the  inviscid  analysis,  equations  (12.52),  to  determine  I 
and  restricting  the  value  of  n  so  that  n  <  2  we  find 


_x 

X 


^  =  r (2-4n)  (|-)1 

o  L  0  J 


1 

2-4n 


(12.57) 


and  the  energy  and  momentum  equations  can  be  written  as  follows 


1  d^e  ^  .  de  _  ,  de  ,  i,  ,,  2  ,«,,2 

R  drj  ' 


1/2  5k  T 


/ 5n-l\ 

\  y  /v-'v 


(^f 


M. 

dV 


(12.58) 


and 


+  ff 


It 


(ifer)  [« - 


=  0 


(12.59) 


Accordingly,  it  can  be  seen  that  since  Moo  =  Moo  (?)  that  the  first  two  terms 

on  the  right  hand  side  of  equation  (12.58)  are  functions  of  I  and  prevent  a  similar 

solution.  If,  however,  the  wall  temperature  is  constant  then  0  ^  ^  (0> 

2 

first  of  these  terms  vanishes.  Finally,  when  Moo  =0  the  second  term  can  be 

2 

neglected  and  a  similar  solution  is  feasible  if  Moo  ^  0  it  can  be  taken  to  be  a 
constant  and  the  similarity  is  then  in  exact.  In  the  present  example  one  can 
clearly  see  that  when  the  flow  is  compressible  the  specification  that  UooOix^ 
is  not  sufficient  to  ensure  a  similar  solution,  but  that  additional  conditions  and 
assumptions  are  necessary. 

Calculations  based  on  the  above  equations  have  been  carried  out  by  Kerrobrock^®^ 
for  Helium  seeded  with  Cesium.  The  free  stream  temperature  and  Mach  number 
were  taken  to  be  3000°K  and  unity  respectively.  A  wall  temperature  of  1500®K 
was  also  chosen. 

Within  the  boundary  layer  being  considered  heat  is  generated  by  two  mechanisms; 
viscous  dissipation,  and  Joule  heating.  Energy  is  transported  from  these  two 
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sources  and  from  the  high  temperature  free  stream  to  the  wall  by  both  conduction 
and  diffusion  of  electrons  against  a  temperature  gradient.  For  the  selection  of 
parameters  made  by  Kerrebrock  both  viscous  dissipation  and  Joule  heating  are  of 
comparable  importance,  and  energy  transport  by  electrons  is  only  a  few  percent 
of  the  overall  heat  transfer. 

To  see  more  clearly  the  influence  of  the  above  heat  sources  on  the  boundary 
layer,  velocity  and  temperature  profiles  are  shown  in  Figures  12-11. 


Figure  12-11.  Temperature  and  Velocity  Profiles  for  Electrode 
Boundary  Layer.  Moo=l,  0^j  =  O.5 


Increasing  values  of  n  correspond  to  free  streams  with  increasing  acceleration. 
This  is  seen  to  correspond  to  the  fact  that  the  velocity  profiles  tend  to  become 
fuller  as  n  increases.  The  temperature  profile  marked  "flat  plate"  corresponds 


to  the  case  for  which  there  is  no  current  flowing  at  all.  The  large  difference  in 

30 

■5—  I  between  the  dashed  curve  and  the  others  is  indicative  of  the  increased  heat 
at]  o 

flux  due  to  Joule  heating  in  the  low  temperature  region  near  the  wall.  For  the 


numerical  example  chosen  the  increase  is  approximately  an  order  of  magnitude. 
The  temperature  excess  at  the  highest  accelerations  can  be  attributed  to  heat 


generated  by  viscous  dissipation  which  cannot  readily  transfer  to  the  wall  due  to 
the  large  amount  of  heat  liberated  near  the  wall  by  Joule  heating. 

The  final  item  of  interest  in  regard  to  the  electrode  boundary  layer  is  the 


electrode  potential  drop.  Due  to  the  fact  that  the  gas  in  the  vicinity  of  the  wall 
is  at  low  temperature  and  consequently  low  electrical  conductivity,  the  electrical 


field  in  the  vicinity  of  the  wall  will  be  much  larger  than  that  in  the  free  stream 


in  order  to  be  consistent  with  a  constant  current  flow.  This  will  lead  to  a 
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larger  potential  drop  across  the  boundary  layer  thickness  than  across  a  corres¬ 
ponding  thickness  of  the  free  stream.  Qualitatively,  the  potential  distribution, 
based  on  the  electric  field  seen  by  a  stationary  observer,  across  a  complete 
channel  will  be  as  shown  in  Figure  12-12. 


Figure  12-12.  Potential  Distribution  Across  a  Channel  of  Height  h. 

The  precise  magnitude  of  the  potential  drop  across  the  boundary  layer  can 
be  calculated  readily.  Consider  the  potential  drop  in  excess  of  the  potential 
drop  through  an  equivalent  thickness  of  free  stream.  Then 


00 


60  =  / 
o 


(E  -  Eoo)  dy 


(Too 


-  B  (Uoo  -  u) 


] 


dy 


(12.60) 


where  it  should  be  noted  that  the  reduction  in  cr  in  the  boundary  layer  will  tend 
to  increase  66  while  the  reduction  in  flow  velocity  there  will  tend  to  reduce  it. 
Depending  on  the  assumed  conditions  60  can  be  positive  or  negative.  Calculations 
of  66  based  on  the  example  cited  earlier  are  shown  in  Figure  12-13. 


n 

Figure  12-13.  Boundary  Layer  Potential  Excess 
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It  is  interesting  to  observe  that  60  is  positive,  for  this  case,  so  that  the 
potential  drop  through  the  boundary  layer  is  indeed  greater  than  in  the  free 
stream  over  the  same  distance. 

Finally,  mention  should  be  made  of  the  magnetohydrodynamic  compressible 
boundary  layer  in  the  region  of  the  stagnation  point  on  a  blunt  body.  The  practical 
example  is,  of  course,  the  re-entering  nose  cone.  For  this  problem  the  inviscid 
flow  has  been  calculated  and  it  has  been  shown  that  Uoo  =  ax  when  is  a  constant, 
the  magnitude  of  a  being  reduced  as  the  strength  of  the  magnetic  field  is  increased. 
These  are  precisely  the  minimum  requirements  for  a  similar  solution,  as  noted 
earlier,  and  with  several  additional  simplifying  assumptions  similarity  solutions 
can  indeed  be  obtained. '  ' 

12.6  MAGNETIC  BOUNDARY  LAYERS 


Up  to  this  point  in  the  present  chapter,  and  in  fact  throughout  this  book, 
problems  of  external  flow  (flow  over  closed  bodies  such  as  airfoils)  have  not 
been  considered.  It  will  be  of  interest,  however,  to  consider  qualitatively  some 
new  boundary  layer  phenomena  that  arise  in  such  flows  when  R^  is  large  and 


R 


is  small. 


(8) 


m 


The  particular  overall  problem  that  will  be  investigated  will  be  the  so  called 

"alligned  flows".  These  are  flows  in  which  the  flow  velocity  and  magnetic  field 

vectors  far  from  the  body  are  parallel.  Now  when  R^^  =  «  and  the  electric  field, 

E,  is  zero  the  Ohm's  law  requires  that  v^B  =  0  so  that  v  and  B  are  not  only 

parallel  at  infinity  but  are  parallel  everywhere.  For  a  body  of  finite  conductivity 

E  =  0  implies  =  0  so  that  the  magnetic  field  is  harmonic  within  the  body. 

Since  it  must  be  a  constant  on  the  surface  the  mean  value  theorem  of  potential 

theory  tells  us  that  it  must  be  zero  everywhere  within  the  body.  Accordingly, 

the  tangential  component  of  B  at  the  surface  must  jump  from  a  finite  value  to 

zero.  As  was  shown  in  Chapter  2  this  corresponds  to  a  current  sheet  at  the 

surface.  Since  v  must  also  go  to  zero  at  the  surface  there  must  also  be  a 

vortex  sheet  there.  When  R^  is  large  but  no  longer  infinite  these  current  and 

vortex  sheets  are  in  reality  boundary  layers.  The  physical  nature  of  such  layers, 

and  the  equations  governing  them  will  be  the  subject  of  this  section. 

®m 


Since,  as  was  seen  earlier,  Pj^  =  the  assumption  of  large  R^  and 
small  Pp  is  tantamount  to  assuming  R  large  and  R  much  larger.  Accordingly 

Ill  “ 

it  can  be  anticipated  that  the  boundary  layer  in  question  will  really  be  two  layers. 

One  will  be  an  outer  layer  in  which  viscosity  is  negligible,  and  the  other  will  be 

-l/2 

a  viscous  sublayer.  The  outer  layer  thickness  will  be  the  order  of  R^  , 
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while  the  inner  layer  will  be  the  order  of  Rg  .  The  flow  external  to  the 
outer  layer,  where  v  aB,  will  be  irrotational.  The  above  physical  interpretation 
of  the  boundary  layer  region  is  illustrated  in  Figure  12-14. 


Figure  12-14.  Sketch  Showing  Inviscid  Magnetohydrodynamic 
Boundary  Layer  and  Viscous  Sublayer 

Before  considering  the  equations  and  specific  boundary  conditions  in  detail, 
it  will  be  of  value  to  first  discuss  the  procedure  for  solution.  First,  the 
potential  flow  must  be  determined  neglecting  both  boundary  layers.  The  boundary 
conditions  for  such  a  solution  are  that  the  surface  be  both  a  fluid  and  magnetic 
streamline.  From  such  a  solution  one  obtains  values  of  u  and  B^  at  the  wall. 
These,  then,  will  serve  as  outer  boundary  conditions  for  the  inviscid  layer. 

The  inner  conditions  for  this  layer  are  evaluated  by  assuming  the  viscous  layer 
to  be  of  negligible  thickness.  Thus,  one  of  the  inner  conditions  on  the  inviscid 
layer  will  be  v  =  0  while  u  0.  Another  inner  boundary  condition  on  the  magnetic 
field  will  be  needed  to  complete  the  formulation  of  the  inviscid  boundary  layer. 

The  viscous  layer  will  then  use  u  (x,  o)  as  obtained  from  the  inviscid  layer  solu- 
ti(m  as  its  outer  boundary  condition  along  with  v  =  0.  The  inner  boundary 
conditions  on  the  viscous  layer  are  the  conventional  ones  of  u  =  v  =  0. 

So  far  little  has  been  said  about  the  boundary  conditions  on  B.  The 
difficulty  lies  in  the  fact  that  when  Rj„  is  no  longer  infinite  then  B  is  no  longer 
zero  within  the  body  and  B  and  B  at  the  surface  are  not  known  prior  to  the 

^  y 

solution  of  the  problem.  The  one  thing  we  do  know,  however,  is  that  when 

R  “  then  B  and  B  should  both  —  0  within  the  body, 
m  X  y 
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Before  resolving  the  above  question  of  the  boundary  values  of  B  and  B  it 

X  y 

will  be  necessary  to  introduce  some  order  of  magnitude  arguments.  It  will  be 
assumed  that  the  thickness  of  the  inviscid  magnetic  boundary  layer,  6.  ,  is 

— 1/2  ^  -I/2 

0  (Rjjj  ),  and  that  the  thickness  of  the  viscous  sublayer,  6^,  is  6  (Rg  '  ). 

Also,  differentiation  with  respect  to  x  will  not  alter  the  order  of  magnitude  of  a 
quantity,  while  differentiation  with  respect  to  y  will  change  the  order  of  magni¬ 
tude  by  6”^.  Whether  it  is  6.  or  6^  will  depend  on  which  layer  is  being  discussed. 

It  was  noted  earlier  that  B„  and  B  tend  to  zero  as  R„  -*  so  that  it  will 

-1/2 

be  valid  to  assume  them  both  6  (Rjj^  )  at  the  wall.  Since  the  viscous  sublayer 
is  extremely  thin  neither  should  vary  appreciably  from  its  value  at  the  wall 
within  this  region.  Next,  consider  the  following  relation 

VxB  =  R^  vxB  (12.61) 

where  v  and  B  have  been  made  dimensionless  by  reference  values  of  v  and  B, 
where  Ohm's  law  has  been  used,  and  E  =  0.  Or,  for  the  present  two  dimen¬ 
sional  problem; 


9B, 


9B 


=  R_ 


9y  ’m 


[uBj-vbJ 


(12.62) 


Within  the  viscous  layer  the  first  two  terms  have  the  following  order  of 
magnitudes 


^  -  S  (R 
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9B 
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9y 


*  =  e 


Also,  within  the  viscous  layer  v  «  u  and  B^  and 


9B 

so  that  9B  /9x  « 

y  9y 

B  are  of  comparable  order  of  magnitudes,  so  that  vB  «uB  .  Accordingly 
X  ^  y 

equation  (12.61)  can  be  simplified,  within  the  viscous  layer  to 

aB„ 


9y 


m  y 


(12.63) 


Or,  integrating  over  the  viscous  layer 


-1/2 


1/2, 


Accordingly,  one  cfui  assume  ®  throughout  the  viscous  sublayer  and  take 

this  as  the  other  inner  boundary  condition  for  the  inviscid  layer. 
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The  equations  governing  both  boundary  layers  will  be  deduced  next,  by 

similar  order  of  magnitude  arguments.  Within  the  inviscid  layer 

aB  aB 

=  0,  so  that  assuming  B^  =  0(1)  leads  to  B^  =  0 

Next,  the  two  momentum  equations  (equation  8.17),  and  equation  (12.62)  are 


considered 


0(1) 

0(1) 

0(1) 

0(1) 

0(1) 

3u 

.  9u 

^  9P  „  ! 

aB 

X 

x  ax  y 

9y 

(12.64) 


av  ,  9v  .  ap  -T 
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ax 


+  B 
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(12.65) 


0(Rm’^^^)  ® 
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uB. 


V  B. 


(12.66) 


B  ^ 

where  P  =  p  +  —5-  and  N  =  °° — 2  R  will  be  assumed  that  P  and  N  are 

Z  jX  p  Uoo 

0(1).  The  resulting  boundary  layer  equations  in  dimensional  form  for  the 
inviscid  layer  are 


au 

H 

+  pv 

9y 

ax 

=  ®x 

9Bx 

ax 

aB 

^  ®y 

(12.67) 

aB 

X 

9y 

=  (TPe 

<uB^ 

-  vBx) 

(12.68) 

ax 

av 
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0 

(12.69) 

aB 

aB 

+ 

ax 

9y 

=  0 

(12.70) 
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with  the  boundary  conditions 

at  y  =  « 

U  =  Uoo  (X) 

* 

®x  ■  ®x«, 

at  y  =  0 

V  =  0 

Bx  =  0 

It  is  interesting  to  see  that  a  boundary  condition  on  at  y  =  0  is  not  needed 

since  when  is  eliminated  from  the  equations  by  equation  (12.70)  they  become 

second  order  in  B  and  there  are  two  boundary  conditions  on  B  available.  The 

value  of  By  at  y  =  0  is  part  of  the  solution  being  sought  and,  accordingly,  it  is 

not  a  suitfjjle  boundary  condition.  From  equation  (12.65)  it  can  also  be  seen 
Q'P  — 1/2  ftp 

that  =  0(Rjj^  )  and  P  =  P(x)  alone.  Accordingly,  can  be  evaluated 

from  the  free  stream  solution. 

The  boundary  layer  equations  within  the  viscous  layer  are  precisely  those 
obtained  earlier  in  the  present  chapter.  The  applied  magnetic  field,  B^.,  is 
taken  to  be  the  value  found  at  y  =  0  from  the  inviscid  boundary  layer  solution. 
The  pressure,  p,  to  be  used  is  P(x)  obtained  from  the  potential  flow  solution. 

hi  conclusion  then,  when  —  «  boundary  layer  phenomena  very  similar 
to  conventional  boundary  layers  can  exist.  They  do,  however,  offer  a  number 
of  new  features  for  future  studies. 
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A  review  of  those  new  phenomena  which  arise  in 
magnetohydrodynamic  boundary  layers  is  presented  in 
[this  report.  The  first  topic  discussed  due  to  its  simplici¬ 
ty  and  relation  to  boundary  layer  flows  is  the  Rayleigh 
problem.  Here  it  is  shown  that  the  introduction  of 
|magnetohydrodynamic  forces  leads  to  basic  changes  in 
the  nature  of  the  flow.  Next  the  basic  boundary  layer 
equations  are  deduced  and  the  conditions  necessary  to 
yield  similar  solutions  deduced.  In  succeeding  sections 
incompressible  and  compressible  boundary  layers  are 
treated  with  attention  paid  to  the  appropriate  external 
boundary  conditions.  Finally,  boundary  layer  phenomena 
arising  when  Rjy^“*  ®  are  considered  in  order  to  illustrate 
some  new  phenomena  within  the  framework  of  boundary 
layer  theory.  It  is  intended  that  the  present  material  will 
be  one  chapter  of  a  forthcoming  book. 


AUTNOft 


cvMnf  IM*  fBCliulB  mm4  ItlAlm  •!!  tli«  tBUltr  Nn*.  Hm 

iVfW  W  Wra  VIVVb 

(JIajC^Cuul  _ _ 


COWNTMIMNW 


